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Abstract. We extend previous work on injectivity in chemical reaction networks to general 
interaction networks. Matrix- and graph-theoretic conditions for injectivity of these systems are 
presented. A particular signed, directed, labelled, bipartite multigraph, termed the "DSR graph" , is 
shown to be a useful representation of an interaction network when discussing questions of injectivity. 
A graph-theoretic condition, developed previously in the context of chemical reaction networks, is 
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applicability of the theory developed. 
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1. Introduction 

Dynamical systems involving networks of interacting elements arise in many fields, 
including chemistry, systems biology, ecosystem modelling, and even beyond the nat- 
ural sciences. The description of these systems requires enumeration of the species 
involved and the interactions between them. Associated with a species is a "state" , 
usually a population or a concentration, and associated with an interaction is a "rate" 
which describes the frequency of the interaction. An important theoretical question 
is: "what claims can be made on the basis only of qualitative knowledge of the inter- 
actions, i.e. without detailed knowledge of functional forms and parameter values?" . 
In this paper we are concerned with questions of injectivity and hence the ability of 
systems to have multiple equilibria. In particular: 

1. Given a wide class of systems (to be defined below) we show how to associate 
with each such system an object termed a DSR graph; 

2. We then show how a simple computation, or even observation, on this DSR 
graph can suffice to tell us that a system cannot have multiple equilibria 
(because the associated vector field is injective). 

The results presented here largely build on those in [U [2] . However the graph- 
theoretic work in this paper is closely connected to two previous strands of work: 
discussions of injectivity /multiple equilibria via interaction graphs [3l|4l[5l[6], and 
discussions of injectivity in chemical reaction networks (CRNs) with reference to so- 
called SR graphs [5] . The classes of systems treated in these references are special 
cases of interaction networks, as defined generally here. 

There is also other important work in this tradition. For example, graph-theoretic 
approaches to questions of monotonicity appear in ^ ilOj and several results in [TT] 
have easy graph-theoretic interpretations. Important questions of persistence in CRNs 
are treated graph-theoretically in [121 US] • 
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The main result is that given any interaction network it is possible to construct a 
multigraph called the directed SR graph (DSR graph) , and check a condition on this 
graph which will rule out the possibility of multiple equilibria. At a practical level, the 
DSR graph is an intuitively meaningful object, closely related to the diagrams drawn 
by researchers in biochemistry and chemical engineering. The condition sufficient for 
injectivity, here termed Condition (*), is simple to state and has been previously 
presented [H [2] . What is remarkable is that a condition originally developed for CRNs 
with mass action kinetics [7], and then generalised to CRNs with looser kinetics [2], 
can actually be applied to arbitrary networks of interacting elements. 

Before proving any results, we will discuss the notion of an "interaction network" , 
both informally and formally. We will see that essentially arbitrary dynamical systems 
on appropriate subsets of M" can be cast as interaction networks, and thus the results 
in this paper have wide applicability. 

1.1. Species and interactions: an informal discussion 

A species S may participate in an interaction R in the following three ways: 

1. Two-way {S ^ R): the species influences the interaction and is itself influ- 
enced by the interaction. 

2. Species-to-interaction (5 — > i?): the species influences the interaction, but 
is unaffected by the interaction. 

3. Interaction-to-species (5 ^ R): the species is influenced by the interac- 
tion, but does not affect the interaction. 

Case (1) describes perhaps the most common situation. During predation, for exam- 
ple, we expect the population of prey both to affect the rate of predation and to be 
affected by it. In case (2) the species is a modulator of the interaction, and this 
modulation may be in a defined direction ("activation" or "inhibition") or an unde- 
fined direction. For example, a simplification common in CRN modelling is to treat 
enzymes as modulators of reactions. In case (3), we say that the species participates 
irreversibly in the interaction. Note that this is slightly different to the usual defini- 
tion of irreversibility for a chemical reaction, because by this definition, if a species S 
is the product of an irreversible reaction and it affects the rate of that reaction, then 
we say that S participates in that reaction reversibly, even though the reaction could 
never run backwards (see for example the model of the TCA cycle presented later). 

A key feature of the DSR graph associated with an interaction network is that 
it has two types of vertices, one representing the species, and another representing 
the interactions. A species vertex S and interaction vertex R may be linked by up 
to two edges, which may be directed or undirected. These encode information about 
how species S affects the rate of interaction i?, and conversely how the dynamics of 
S are affected by the rate of interaction. A formal treatment will be developed later, 
but the reader might like to glance at Appendix [Cl to get a feel for the meanings of 
various motifs in DSR graphs. 

1.2. Summary of the results 

Rectangular domains. A subset X of M" is a rectangular domain iff for i = 
1, . . . ,n there exist quantities — oo < Xi^min < oo and —oo < Xi^max < oo satisfying 

Xi^min ^ Xi^maxj and for Cacll X = [xi, . . . , Xn\^ G X 

^i,min ^ •^i ^ •^i,max ^ 1, . . . , 77 . 

Here the symbol -< can mean < or <. For example X could be all of R", an orthant, 
etc. Note that X may be of lower dimension than M". 
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Injectivity and multiple equilibria. A function f : X ^ Y between arbitrary 
spaces X and Y is injective iff whenever a,b G X , a b, then /(a) ^ /(&)• Let 
X C M", and f : X M" an arbitrary (i.e. continuously difFerentiable) function 
on X. Supposing some model of the natural world gives rise to a vector field x = /(x), 
or the map Xn+i = Xn-\- f{xn), where x G X. Suppose further that we have incomplete 
data, so that all we know is that / belongs to some set of functions J-^ but by some 
means we can show that aX\ f G T are injective. This immediately implies that each 
f G T can have no more than one zero, i.e. the vector field can have no more than 
one equilibrium, or the map no more than one fixed point. 

The key idea: decomposing functions. Our central theme in this paper is the 
following: we start with a set of functions from some rectangular domain X C M" 
to R". Now for some integer m > 1, and for each f G T we write f = f2 ° fi, where 
fi : X R™ and /2 : M™ R". Such decompositions always exist (for example 
we can choose m = n, fi = f, and /2 to be the identity, or vice versa), but often 
there are natural decompositions implied by the way that the models are constructed. 
Intuitively, we can think of X as "species space" , the set of allowed species values, 
and R™ as "interaction space" , the space of allowed interaction rates. Differentiating 
/ = /2 ° /i gives Df{x) = Df2{fi{x))Dfi{x), i.e. the Jacobian Df has a natural 
product structure. Specifically, {Df2)ij encodes information about how interaction 
j affects affects quantity i, while (Dfi)ji encodes information about how quantity i 
affects the rate of interaction j. 

Pq matrices. The product structure in the Jacobian forms the theoretical 
starting point for the subsequent treatment. In particular, the class of -Pq "* matrices 
(defined in Appendix |^ is shown to be a useful class to consider when asking ques- 
tions about injectivity. Two very general lemmas are developed: Lemma 13.31 gives a 
sufficient condition for any set of matrices with product structure to be -Pq""* matrices; 
Lemma 13.41 provides restrictions under which the sufficient conditions in Lemma 13.31 
are also necessary. 

DSR graph results. While the matrix-theoretic approaches give sharp results, 
there are very elegant and intuitive graph-theoretic results which arise as corollaries 
of these. Roughly speaking, the DSR graph is constructed as follows: returning to 
the decomposition of the Jacobian Df{x) = Df2{fi{x))Dfi{x), the elements {Df2)ij 
and {Dfi)ji together contain information on how species i is affected by, and affects, 
interaction j, and unless both are identically zero, together they map to one or more 
edges between species vertex i and interaction vertex j in the DSR graph. Via a series 
of lemmas we find that a previously developed condition on cycles in the SR graph, 
restated for the DSR graph, guarantees injectivity of the interaction network with 
some outflow/degradation of each species. The implications for multiple equilibria 
under different assumptions are summarised in Corollary 14. 2| which can be seen as 
the main result in this paper. 

1.3. An example: the TCA cycle 

Before presenting the theoretical development, or even rigorously defining interac- 
tion networks and DSR graphs, it is helpful to present a nontrivial example. Further 
discussion of this and other examples will be provided after development of the theory. 

We consider the model of the TCA cycle discussed in [14 . The model is a fairly 
complicated biological model, and is presented in stages. At each stage the relevant 
question is whether it can admit "multiple positive nondegenerate equilibria" , a notion 
to be made precise later. 

The backbone of the model is a cycle of eight interconversions, that is, very 
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basic chemical reactions in which one substrate is simply converted into anothei0. 
Importantly, some of the reactions are assumed to be irreversible while others are 
reversible. The DSR graph for the model at four different levels of complexity is 
shown in Figure [TTT] 




Fig. 1.1. The DSR graph for the model of the TCA cycle presented in \14^ at four different 
stages of construction. Details of the biochemistry can be found in Negative edges are repre- 

sented with dashed lines, while positive edges are represented with bold lines. Apart from the one 
edge labelled oo all edges have edge-label 1. Various quantities which do not affect the cycle structure 
have been omitted from the DSR graphs. Implications are described in the text. 



1. The basic structure of the model gives the DSR graph shown in Figure fTTTb . 
Models with this simple structure cannot display multiple positive nonde- 
generate equilibria. In fact previous theory [15] indicates that with mild 
additional assumptions, there must be a unique, globally stable, equilibrium. 

2. Adding the inhibition by oxaloacetate of succinate dehydrogenase gives rise 
to the DSR graph shown in Figm-e fTTTb . Analysis of this graph tells us that 



^The analysis remains the same if we treat it as a cycle of nine interconversions, including cis- 
aconitate. 



the conclusion about multiple nondegenerate equilibria remains true. 

3. Several reactions in the TCA cycle cause reduction of NAD. As [NAD] + 
[NADH] is a conserved quantity, one of the pair can be eliminated. Including 
NADH gives the DSR graph shown in Figure fTTTb . Again, the model does not 
allow multiple positive nondegenerate equilibria, despite the large number of 
cycles in the DSR graph. 

4. Finally, adding the aspartate amino transferase (AAT)-catalysed reaction ef- 
fectively adds an extra interconversion between a-kctoglutarate and oxaloac- 
etate giving rise to the DSR graph shown in Figure fTTTH . Now the DSR graph 
contains "bad cycles" and it is no longer possible to rule out multiple positive 
nondegenerate equilibria from the graph. 

2. Interaction networks 

Assume that there are n species with "amounts" (concentrations, populations, 
etc.) xi, . . . ,Xn, and define x = [xi, . . . ,a;„]"^. Let there be m interactions which 
occur with rates vi{x), . . . , Vmix), each involving any subset of the species, and define 
v{x) = [vi{x), . . . ,Vm{x)]'^ ■ Finally, define the ith interaction function, fi{v{x)), 
to be the total rate of production/consumption of species i. The evolution of the 
system is then given by: 

±^^f^{v{x)) l = l,...,n. (2.1) 

It is assumed that all functions are and that each Xi takes values in some interval 
(bounded or unbounded), so that the state space is a rectangular domain in R". Two 
absolutely general features of (|2.1[) are: 

• By the chain rule, the Jacobian (at each point in space) allows a decomposi- 
tion J — SV where Sij — and Vji = -g^. 

• The interaction structure can be represented as a bipartite multigraph (dis- 
cussed informally above, and to be defined formally later) with n species 
vertices and m interaction vertices. 

A slight, but important, variant on Equation (j2.ip involves assuming nonzero 
outflow or degradation rates for each species. We introduce a set of scalar functions 
qi(xi), i = 1, . . . ,n, each satisfying > throughout its domain of definition. The 
system becomes 

ii = fiiv{x)) - qi{x^) i^l,...,n. (2.2) 

The Jacobian is now J — SV — D where D is the positive diagonal matrix defined by 
Du = ^^. Defining 

f{v{x)) = [/i(u(a;)),...,/„(u(x))]^ and Q{x) = [gi (xi ),..., g„(2;„)]^ , 

we can write (|2.1I) and (|2.2p in abbreviated forms which we term \Nq\ and \N+\ : 

X = j(v{x)) (iVo) 

and 

x = j{v{x))-Q{x). {N+) 

The situation where some components of Q are identically zero is also important. 
Given any 9 C {I, . . . ,n} we can define a class of systems which we term jNel : 



X = f{v{x)) - Qe{x), 
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with Qg{x) — [qi{xi), . . . , qn{xn)V , <li{xi) — for i ^ 9 and ^ > for i G 0. Note 



that both dA^oD and (N^ ) are now special cases of (j^J), with 6 — 9 and 9 — {1, ... ,n} 
respectively. 

2.1. Arbitrary dynamical systems as interaction networks 

Consider a rectangular domain X C K", and some function F : X ^ M". 
With Q{x) and Qg{x) defined as above, consider the dynamical systems: 

x^F{x), x^F{x)-Q{x) and x = F{x) - Qe{x). (2.3) 

Defining id(-) to be the identity on R", we can certainly write = F o id, and the 
three systems can be rewritten 

i = F(id(x)), X ^ F{iA{x)) ~Q{x) and x ^ F{id{x)) - Qg{x). (2.4) 



Thus the three systems in (|2.3p can be cast as instances of \N+\ and \Ng\i , and 

all theory developed in this paper can be applied. Of course the choice = o id is 
not unique. We will see by example later that this formal approach may be useful, 
but does not necessarily give the strongest results. This theme will be more fully 
developed in future work. 

2.2. Linear interaction functions 

If the functions /,; take the form 

m 

fz{v{x)) = ^S^jVj{x), 



for some constants Sij, then \Na\ , \Nj^\ and \Ng\ reduce, respectively, to 

X = Sv{x), X — Sv{x) — Q{x) and x = Sv{x) — Qg{x). 

These equations are familiar: they are the equations (with or without outflow) for 
a general CRN. In this context, x becomes the vector of reactant concentrations, 
V is the vector of reaction rates, and S is the "stoichiometric matrix". However 
the assumption of linearity is also common beyond CRNs, for example in models of 
regulatory networks (such as that in [16 , discussed in the examples later). 

3. Matrix-theoretic results 

3.1. Sets of real numbers and matrices 



This paper is concerned with what can be said about the dynamics of < \No\i , ( 
and IjNgl , knowing only that S belongs to some matrix-set S, and V belongs to some 
matrix-set V, possibly related to S. For this reason we start by developing notation 
and ideas connected with sets of matrices. Definitions and notation which are well 
known are summarised in Appendix [Xl 

Define M>o = (0, oo), R<o = (-oo, 0), M>o = [0, oo) and R<q = (-oo, 0]. A set of 
real numbers TZ is signed HTZ C K>o or 72. C M<o or TZ — {0}, and is weakly signed 
if 7?, C M>o or 72 C M<o. A set of real numbers which fails to be weakly signed (i.e. 
intersects both K>o and M<o) is unsigned. 

When A4 is some set of matrices or real numbers, the statement "Af has property 
P for all M G A^" will be abbreviated to "Ai has property P". So, for example, the 
statement "det(A() = 0" should be read as "det(M) for each AI G A(". 
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Sums and products of matrix-sets. Given two sets of matrices, A and B, we 

will always assume that A and B arc defined as the ranges of matrix- valued functions 
A, B over some underlying space X, so that A = {A{x)\x G X}, B = {B{x)\x £ X}. 
This entails no loss of generality, because even where the matrix-sets are defined 
purely set-theoretically, or via functions over different spaces, it is an easy matter to 
redefine the sets as the ranges of matrix-valued functions over some common space. 
Suppose, for example, A is defined as the set of all n x n matrices, and T> as the set 
of all n X n positive diagonal matrices. In this case we could define X = R"^" x IR"g, 
with A now the range of a function A taking the first coordinates of X to elements 
of A, and V the range of a function V taking the final n coordinates to the diagonal 
elements of "D. 

Now as a convenient abbreviation we define AB — {A{x)B{x)\x G X} and A + 
B = {A{x) + B{x)\x G X}. If, for each A & A,B € B, there exists x € X such 
that A = A{x),B = B{x), then we will say that A and B are independent. For 
independent sets of matrices, it follows that AB = {AB\A G A, B G B}, and similarly 
A + B = {A + B\A G A,B G B}. Thus if two sets are independent then their 
sum/product is their pointwise sum/product. 

As an example, let x = [x-i,X2,xs,X4]'^ G M'* and for each x define 



A(x) 



X\ X2 
X3 X4 



and B{x) 



xi 2 

-X3 1 



X2 
X4 



Here A = {A{x)\x G M^} and B = {j3{x)\x G M^} are clearly not independent. In 
fact, A + B = {A{x) -h 13{x)\x G M^} consists of the single matrix 



1 2 
1 



On the other hand, suppose 



A{x) 



Xi 
X2 



and B{x) = [ X3 — xi X4 — X2 ] 



In this case A = {A{x)\x G M.^} and B = {B{x)\x G M*} are independent, and it is 
easy to confirm that AB = {AB\A & A,B & B}. 

Notation. Let M be an n x m matrix, with 6 C {1, . . . , n} and 7 C {1, . . . , to}. 
The following notation will be used: 

• M{S\j) is the submatrix of M with rows indexed by S and columns indexed 
by 7. 

• M[5|7] = det(M(^|7)). We write M[5] as shorthand for M[5|(5]. 

• Ms-y is an n X TO matrix defined by {Ms^)ij — Mij if i <E S and j G 7 and 
{Msj)ij = otherwise. Note that all \5\ x I7I submatrices of M^^, apart 
possibly from M((5|7), must contain a row or column of zeros (and hence, if 
they are square, must be identically singular). 

A set of real n x to matrices can be seen as a subset of M"^™ and thus inherits 
topological properties such as openness, connectedness, etc. The closure of a matrix- 
set is denoted by c\{A4). and / denote the zero and identity matrices with 
dimensions being clear from the context. 

Entries in matrix-sets and minors of matrix-sets. Given a matrix-set 
A = {A{x)\x G X}, we define Aij : X ^ R by Aij{x) = {A{x))ij, and define Aij = 



{Aij{x) I X e X}. A product of entries such as AijAki means {Aij{x)Aki{x) \ x e X}, 
and this notation extends to arbitrary products and sums of products of entries. 
Given matrix-sets A = {A{x)\x G X} and B = {B{x)\x e X}, A[S\^]B[j\S] should be 
interpreted as {.A(a;)[^|7];B(x)[7|(5] \ x S X}. Expressions involving sums/products of 
entries/minors always follow these conventions, which considerably simplify notation. 

Sign-classes. A set of n x m matrices Ai will be defined as a sign-class if given 
any M G M, and any S C {1, . . . , n}, 7 C {1, . . . , m}, Ms-y G c\{M). Since 7 and 6 
may both be empty, cl(A^) contains the zero matrix. 

Sign-classes may have related elements. For example, let x = [xi, X2, X3, X4]'^ G 
and define 



Mix) 



X1X3 X1X4, 
X2X3 X2X4 



M = {M{x)\x eE."^}. 



It is easy to check that M is a. sign-class, even though the entries are not independent, 
and there is no M G d{M) with Mn = but M12, M21 ^ 0. 

Determinant expansions. Given a set of n x n matrices Ai = {Ai{x)\x G X}, 
and a permutation a of the ordered set [1, . . . define P{a) to be the parity of 
the permutation, i.e. P{a) = 1 for an even permutation and P{a) = — 1 for an odd 
permutation. Corresponding to a there is a term in the determinant expansion of 
M, Ta. Define f„ : X ^ R, by fa{x) = P{a) HlLi (Ma;))*^. , and = {fo,{x)\x G 
X}. The conventions adopted for sums and products of matrix-sets apply to terms in 
determinant expansions as well, so it is possible for two terms to sum to zero, even 
though neither consists of a singleton. As an example, let x = [xi,X2]'^ G and 
define 



Mix) 



Xi X2 
Xi X2 



M = {M{x)\x G M^}. 



Then for a the identity and /3 — (21), — {2:1X2 | Xi, a;2 G M} = M, T/j = 
{—X1X2 \xi,X2 G M} = M; however + = {a;ia;2 — X1X2 \xi,X2 = {0}. 

3.2. systems 

Our concern is to find conditions which tell us when an interaction network for- 
bids multiple equilibria. This is the case if the system is injective, i.e., the right hand 
side of the dynamical system is injective. We cannot in general expect jNgl to be 
injective: quite generally we may get conserved quantities and hence invariant mani- 
folds foliating the state space (termed stoichiometric classes for CRNs). Often these 
manifolds are affine subspaces, that is, they are cosets of some linear subspace. Define 
the conditions: 

CO. ( pVol ) have P^^^^ Jacobians. 

have Jacobians. 
have nonsingular Jacobians. 
is injective on any rectangular domain. 



A, 



A, 



CI 
C2 
C3 
C4 

any invariant afhnc subset of the phase space. 



Ng I has no more than one nondegenerate equilibrium in the relative interior of 



By basic results on Pq matrices, CO implies CI and hence, trivially, C2. Below, 
in CoroUarv 13.21 we show that C2 implies CO, and so CO, CI and C2 are equiv- 
alent. Techniques such as those employed in [51[T71[I3 to prove absence of multiple 
equilibria using degree theory require only a nonsingular Jacobian rather than P or 
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P^^-' Jacobians. The equivalence of CO, CI and C2 shows that this requirement 
is only apparently weaker for systems which may have arbitrary outflows. It is well 
known that if a function / has a P matrix (or p(~) matrix) Jacobian on some rect- 
angular domain in M", then / is injective on this domain (Theorem 4 and subsequent 
remarks in [20]), a fact which was applied in [4l[T]. Thus CI implies C3. Finally, 
C3 implies C4. The meaning and proof of this fact are presented in Appendix [b1 
Thus CO implies all the other statements. Defining a "Pq ^ system" to be a system 
with Pg Jacobian, we have: 



If dA^oD is a Pq — ) system, then ^N+l is injective, and ( [TVeD admits no 



multiple nondegenerate equilibria on the relative interior of any in- 



variant afflne subset of the phase space. Note that ( A^_|_ ) was obtained 
from ( |7VoD by assuming that each Xi is subject to outflow or degrada- 
tion, and jNel was obtained from (Nq) by assuming that some Xi's 
may be subject to outflow or degradation. 
Remark. Note that in the case of some chemical reaction systems, inectivity of 



( ) implies injectivity of ^Ngj for any 9 

Lemma 3.1. Consider any n x n matrix A and let V be the set of n x n positive 
diagonal matrices. If A fails to be a Pq matrix, then A — V = {A — D \ D £ P} 
contains matrices with determinants of all signs. 

Proof Given any D eV, write A- D = -D{I + D-^{-A)), so that sign(det(A- 
D)) — — sign(det(/ + _D^^(— A))). Expanding, we have 

det{I + D-\-A))^l+ iD-\-A)m 

<5C{l,...,n} 

where the sum is taken over all nonempty subsets of {1, . . . , n}. Using the Cauchy- 
Binet formula [21] , and noting that the only nonzero minors of a diagonal matrix are 
the principal minors, gives {D~^{—A))[S] = D~^[S]{—A)[6]. 

First, for any e > define ^'(e) £ P by D{e)u = 1/e. Then [D{e)-% = e, and 
iDie)-^)[d] =el*l. Sodet(/ + Z?(e)-i(-A)) = l + ^5^{i^ „jel*l(-A)[5]. For small 
enough e > 0, clearly det(/-|-_D(e)^^(— A)) > 0. Thus there exists a positive diagonal 
matrix Di such that sign(det(A — Di)) = —1. 

Suppose that A is not a Pq""* matrix. This means that there is some nonempty 
So C {l,...,n} such that {-A)[So] < 0. Now for any e > define D{e) e 13 by 
D{e)u = e if i e (5o and D{e)^^ = 1/e if i ^ <5o. In this case {D(e)-^)[S] = gl^WI-l^n-Sol^ 
From this we get that 

det{I + D{e)~\-A))^l+ ^ el5\^o|-|5n5ol(_^)[^] 

<5C{l,...,r!} 

= e~^^°^{—A)[So] + higher order terms in e. 

Clearly, for small enough e > 0, det(/ + D{e)^^{—A)) < 0. Thus there exists a 
positive diagonal matrix D2 such that sign(det(A — D2)) = 1. By continuity of the 
determinant, there exists A G (0,1) such that det{A — XDi — {1 ~ X)D2) — 0. As 
ADi + (1 - X)D2 G P, the result is proved. □ 

Corollary 3.2. Consider any set of n x n matrices A and let V be the set of 
n X n positive diagonal matrices. If any A (z A fails to be a Pq ' matrix, then 
A — T) = {A — D\AeA, De T>} contains matrices with determinants of all signs. 
Proof. This is an immediate corollary of Lemma 13.11 □ 
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Remark. Together with the basic fact that a Pq matrix plus a positive diagonal 
niatrix is a P matrix, Lemma [3.11 impHes the following characterisation of Pq matrices: 
a matrix A is a Pq matrix iff det(j4 + D) > for every positive diagonal matrix D. 

3.3. Identifying matrix products as Pq ^ matrices 

Having seen that whether a set of matrices are Pq^^ matrices is central to ques- 
tions concerning multiple equilibria, the next stage is to present a general rule for 
deciding when the product of two matrix-sets iS and V (independent or otherwise), 
consists entirely of Pq matrices. A very general sufficient condition is provided in 
Lemma 13.31 This condition is also proved to be necessary in Lemma |3.4[ provided S 
and V are independent and at least one of 5 or V is a sign-class. 
Lemma 3.3. Consider any set ofnxm matrices S and any set ofmxn matrices V. 
Assume that for any (5 C {1, . . . , n}, 7 C {1, . . . , m} satisfying \j\ = \5\, S[S\j]V['y\S] C 
(— 1)I''IM>o. Then SV consists of Pq matrices. 

Proof. The result follows from the Cauchy-Binct formula. In the usual way, let 
iS = {S{x) I X £ X}, and V — {V{x) \ x G X}. Given some x € X and some S = S{x) 
and V ~ V{x) we have 

{svm^ ^ smv[i\5]- 

7C{l,...,m} 
l7l = |5| 

Supposing the conditions of the lemma are fulfilled, then S'[(5|7]y[7|(5] € (— 1)I''Im>o, 
and thus {SV)[d] G (-1)I'^'Im>o. This proves that SV is a Pq^"' matrix. □ 

Remark. If S and V are independent, the condition iS[(5|7]V[7|(5] C (— 1)I''IR>o 
means that either, i) iS[(5|7] = 0, or ii) V[7|(5] = 0, or iii) iS[i5|7] and V[7|(5] are both 
weakly signed with either S[S\-f], {-V)[-f\6] C M>o or S[S\j], (-V)[7|^] C K<o. 

Depending on S and V, the conditions in Lemma [531 mav be necessary as well as 
sufficient to guarantee that SV consists of Pq ' matrices. In particular the following 
case where either 5 or V (or both) are sign-classes of matrices often arises: 
Lemma 3.4. Consider any set ofnxm matrices S and any set of m x n matrices 
V. Assume that S and V are independent, one of S or V is a sign-class, and SV 
consists of Pq matrices. Then S[S\j]V[j\S] C (— 1)''''M>o for each 5 C {!,..., n}, 
7 C {1, . . . ,m} satisfying \ — \S\. 

Proof Suppose there are some sets S and 7 such that iS[^|7] V[7|(5] ^ (— l)l'^llR>o- 
This means that there exist S eS andV eV such that S'[(5|7]y[7|(5] e (-l)l'^l+iM>o. 
For definiteness assume that V is a sign-class of matrices, so that Vys 6 cl(V), and 

so {SV^sM = S[d\-/]V[j\6] e (-l)l'^l+iK>o. Thus SV^s fails to be a Pq*"' matrix. 
Since the class of Pq ' matrices is closed, its complement is open. Thus any matrix 
sufficiently near to SV-ys fails to be a Pq ' matrix. Since, by independence of S and 

V, SV^s £ cl{SV), there are matrices in SV which fail to be Pq"'- The argument 
works equally well if iS is a sign-class. □ 

There are a couple of immediate observations to be made from the above lemmas. 
Lemma [3.41 places the following basic restriction on systems with unsigned entries and 
Pq ^ Jacobians: 

Corollary 3.5. Consider a system defined by matrix-sets S and V fulfilling the 
conditions in Lemma \3.4\ If for some i 6 {1, . . . ,n}, j £ {!,..., m}, Vji is unsigned, 
then Sij — 0. Similarly if Sij is unsigned, then Vji = 0. 

Proof. Setting S — {i} and 7 = {j} in Lemma [3.41 gives us the result. □ 
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Loosely speaking, this corollary tells us that in a Pq system, if a quantity has 
influence of unknown sign on an interaction, then it should not itself be affected by the 
interaction. In particular, in CRNs, if a substrate occurs on both sides of a reaction 
with unknown influence on the rate, then the substrate must occur with the same 
stoichiometry on both sides. As an example of how more complicated "forbidden 
scenarios" can also be formulated, we have the following result: 

Corollary 3.6. Consider a system defined by matrix-sets S and V fulfilling the 
conditions in Lemma \3.4\ with V a sign-class. Suppose for some indices k, Sji ^ 
and Vkj is unsigned. Then there is no index I such that Sik 7^ and Vu 7^ 0. 
Proof. Assume the contrary and consider the 2x2 submatrices 



From Corollary 13.51 Sjk must be zero, and by assumption Sji, Sik 7^ 0, and so 
A:}] 7^ 0. Further, Vu 7^ 0, and Vkj is unsigned, so, since V is a sign-class, 

V[{i, fc}|{j, I}] is unsigned. By Lemma [3^ SV are not all Pq matrices. □ 

The graphical presentation of this result in Section 16.31 provides some intution as 
to its meaning. 

Remark on previous matrix-theoretic results. The special case described 
in [ij [2] of chemical reactions with no substrate on both sides of any reaction (termed 
NAC systems in the first reference and NIC system in the second) corresponds to 
S consisting of a single matrix, i.e. S = {S}, and V = Qq{—S'^) (see Appendix [X] 
for a definition). In that case the condition iS[(5|7] V[7|<5] C (— 1)I*IR>o in Lemma [331 
implies that all square submatrices of S are either singular or sign nonsingular (see 
Appendix Since {5*} and Qq{—S^) are trivially independent, and Qo{—S'^) is 
a sign-class of matrices, this condition on S is both necessary and sufficient for all 
matrices SV to be Pq •* matrices. In 1^ there was also some discussion of a model 
where S was not a single matrix. 

Lemmas 13.31 and 13.41 directly give a number of easy generalisations of matrix- 
theoretic results in ^ li which will be developed in future work. Here we concentrate 
on the graph-theoretic corollaries. 

4. SR graphs and DSR graphs 

4.1. SR graphs 

SR graphs for single (rectangular) matrices. SR graphs (or species-reaction 
graphs) are bipartite multigraphs with two vertex-sets termed S-vertices and R- 
vertices, and signed, labelled edges. Although originally defined for CRNs '7], they 
can also directly be associated with matrices as in pj: given any rectangular matrix 
M, we can construct the associated SR graph Gm where an edge exists between S- 
vertex i and R-vertex j iff Mij ^ 0. The edge takes the sign of Mij and the edge-label 
\Mij\. When the matrix M is the stoichiometric matrix of a CRN with SR graph G, 
then Gm — G for any reasonable kinetics iff no substrates occur on both sides of any 
reaction in the system. This special case is treated in detail in 2]. 

SR graphs for sets of (rectangular) matrices. Given any matrix-set M., 
we can construct an SR graph Gm which is, roughly speaking, the amalgamation 
of the SR graphs Gm associated with each M £ ^A. If Mij = {0}, then there is 
no edge between S-vertex i and R-vertex j; if A4ij 7^ {0} and Mij C [0, 00), then 
there is an edge between S-vertex i and R-vertex j with sign if Mij 7^ {0} and 
Mij C (— cx),0], then there is an edge between S-vertex i and R-vertex j with sign 



Si{j,l}\{i,k}) 




and V{{i,k}\{j,l}) 



Vij Vu 
Vkj Vki 
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— 1; finally if Mij is unsigned then we introduce a pair of oppositely signed edges 
between S-vertex i and R-vertex j. If A4ij = {k} where k ^ 0, then we give the 
unique edge between S-vertex i and R-vertex j an edge-label of |fc|, otherwise we give 
all edges between these vertices edge-labels of oo. An example of a matrix-set and 
the associated SR graph is shown in Figure [4TT] 



M = 



a X 
2 b 



Si 



■(> 



i?2 



■Ri 



■S. 



Fig. 4.1. Assume that a,b can take all values in [0,1] and x can take all values in [—1,1]. 
The SR graph associated with M contains five edges. Positive edges are represented with bold lines, 
negative edges are represented with dashed lines. These conventions will be followed throughout. 



Cycles in SR graphs. Cycles in SR graphs are minimal undirected paths from 
some vertex to itself. As seen in Figure |4?T1 an SR graph may have a pair of oppositely 
signed edges between S-vertex i and R-vertex j forming a cycle of length 2. Such cycles 
will be termed short cycles. Cycles of length greater than 2 are then long cycles. 
Since all edges in an SR graph are signed, all cycles in an SR graph have a sign, 
defined as the product of signs of edges in the cycle. A cycle C in an SR graph also 
has a parity P{C): it is either an o-cycle or an e-cycle according to whether 

P(C) = (-l)lc^l/2sign(C) 

is negative or positive. Note that short cycles in an SR graph are always e-cycles. 
Given an edge e, define val(e) to be the edge-label of e. When C is a cycle containing 
edges ei, 62, ... , e2r such that and e(i mod 2r)+i are adjacent for each i = 1, . . . , 2r, 
we can define its stoichiometry as follows: if any edge in C has edge-label 00, then we 
set stoich(C) = oo; otherwise 



stoich(C) 



]^val(e2i-i) - ]^val(e2i) 



Note that this definition is independent of the starting point chosen on the cycle. A 
cycle with stoich(C) = is termed an s-cycle. An e-cycle which is also an s-cycle is 
an es-cycle. A disconnecting partition of a cycle C is the (unique) partition of 
the cycle into two edge-sets, such that no two edges in either set share a vertex. 

S-to-R intersection in SR graphs. The intersection of two cycles in an SR 
graph can be divided into a set of vertex-disjoint paths. We say that two cycles have 
S-to-R intersection if each component of their intersection is an S-to-R path, i.e. a 
path between an S-vertex and an R-vertex. 

4.2. DSR graphs 

Where SR graphs are associated with matrix-sets, DSR graphs are associated with 
pairs of matrix-sets. Given any two sets of n x m matrices, A and the DSR graph 
Qa,B is a signed, labelled, bipartite multigraph with edges which may be directed 
or undirected. Here the definition is presented, while a more intuitive discussion is 
presented in Appendix [Cl 
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First, we take the SR graph Gj^ and create a directed version by insisting 
that all edges are directed from R- to S- vertices. Similarly, we define Gg, a directed 
version of Gs with all edges pointing from S- to R-vertices, and all edge-labels set to 
be OQ. ^ _^ 

1. Since S- vertices in and Gg can be identified, and similarly for R-vertices, 
Ga and Gb can be amalgamated into a single sj^ned, labelled, directed 
multigraph G with all edges from both G^ and Gg. (Any S- vertex and 
R-vertex in G may be connected by 0, 1, 2, 3, or 4 edges.) 

2. If two vertices in G arc connected by a pair of edges with the same sign 
but opposite orientation, then we replace these with a single undirected edge 
with sign and edge-label imported from G^. The resulting graph G^^g is 
now termed the DSR graph. Any S-vertex and R-vcrtc;x in G^^g may be 
connected by 0, 1 or 2 edges, some of which may be directed. 

Every edge-set in G^ has a corresponding edge-set in G^^g with the same edge- 
signs and edge-labels. Thus ignoring direction of edges, Gj, is a true subgraph of 
Gj^^B- Similarly, ignoring direction and labels on edges, Gb is a true subgraph of 
Ga,b- Given matrix-sets A and B, Gj^^b is not in general the same as Ge_^, although 
if we ignore edge-labels, reversing the orientation of all directed edges takes G^^g to 
Gb,a- 

Subgraphs of DSR graphs. If a DSR graph G is associated with a pair {A, B) 
of sets of n X m matrices, then given S C {1, . . . , n} and 7 C {1, . . . , m}, G{6\j) will 
refer to the subgraph of G associated with the pair {A{S\'y),B{6\'y)). DSR subgraphs 
with an equal number of S- and R-vertices will be referred to as square DSR graphs. 

Directed paths in a DSR graph. An edge in a DSR graph has R-to-S 
direction if it is undirected or directed from an R-vertex to an S-vertex. It has S- 
to-R direction if it is undirected or directed from an S-vertex to an R-vertex. So 
undirected edges have both R-to-S direction and S-to-R direction. A subset of the 
edge-set of a DSR graph has R-to-S (S-to-R) direction if all edges in the set have 
R-to-S (S-to-R) direction. Two edge-sets in a DSR graph are oppositely directed if 
all edges in one have R-to-S direction and all edges in the other have S-to-R direction 
(one or both sets may contain undirected edges). A directed path in a DSR graph 
consists of alternating S-to-R and R-to-S edges. 

Cycles in DSR graphs. A cycle C in a DSR graph is a minimal directed path 
from some vertex to itself. Of course, some edges in a cycle may be undirected. Cycles 
are always either e-cycles or o-cycles, since all edges are signed. They may also be 
s-cycles. The following two remarks follow immediately from the definition: all short 
cycles in a DSR graph are e-cycles but not es-cycles because at least one edge has 
edge-label 00; any cycle has a disconnecting partition consisting of an S-to-R edge-set 
and an R-to-S edge-set of equal size. 

Formal cycles. Given a DSR graph G we can construct an undirected version 
of it, G by making all directed edges in G undirected in G . Cycles in G will be 
termed formal cycles in G and may or may not be cycles. Sometimes we will refer to 
a cycle in G as a "genuine" cycle to emphasise that the cycle is directed. 

Orientation of cycles. If a (genuine) cycle G in a DSR graph contains only 
undirected edges, then it has two natural orientations. On the other hand, if G 
contains some edge which fails to have both S-to-R and R-to-S direction, then C has 
one natural orientation. Thus we can always choose cither one or two orientations 
for any genuine cycle. If we choose an orientation for a cycle G, then each edge 
(including undirected edges) in G inherits an orientation, which we can call that edge's 

13 



"C-orientation" . We say that two cycles C and D have compatible orientation if 
we can choose an orientation for C and an orientation for D such that each edge 
in their intersection has the same C-orientation and Z3-orientation. Note that two 
cycles with empty edge-intersection trivially have compatible orientation, and that 
two undirected cycles may be incompatibly oriented. 

S-to-R intersection between cycles. Two cycles in a DSR graph are said to 
have S-to-R intersection if they have nonempty intersection, compatible orientation, 
and moreover each component of their intersection has odd length, i.e. it is either an 
S-to-R path or an R-to-S path. (Note that there is no implied direction in the term 
"S-to-R intersection"). 

4.3. The main results 

We have now developed sufficient terminology to state the main results of this 
paper. Define: 

Condition (*): All e-cycles in a DSR graph are s-cycles, and no two e-cycles 
have S-to-R intersection. 

Theorem 4.1. Consider a set of n x m matrices S, a set of m x n matrices V, and 
the corresponding DSR graph G = G^^^t. If G satisfies Condition (*), then S and 
V satisfy the conditions in Lemma \3.3i 

The proof of this theorem will be presented after considerable machinery has been 
built up. The following corollary, which is the main result in this paper, tells us what 
Theorem 14.11 implies. 

Corollary 4.2. Consider an interaction network jNoD defined on a rectangular 
domain X C M", with DSR graph G satisfying Condition (*). Then the Jacobians of 
UNqI ) are all Pq matrices. Associated systems of the form have Jacobians which 

are all P(-) matrices, and hence cannot have multiple equilibria in X. Associated 
systems of the form fjVe[ j, including fiVpl j itself cannot have multiple nondegenerate 
equilibria in the relative interior of any invariant affine subset of X . 

Proof. These claims follow immediately from Theorem 14.11 coupled with the dis- 
cussion in Section [321 D 

Remark. Note that there is an alternative "pointwise" phrasing of Corollarv l4.2l 
Rather than associating a single DSR graph G with the interaction network ( |iVoD , we 
can associate a DSR graph G{x) with each point x ^ X. Provided each G{x) satisfies 
Condition (*), then the Jacobians of IjNo]) are all Pq •* matrices, and all the other 
conclusions follow. It can occur that G fails Condition (*), but at each x & X, the 
DSR graph G{x) satisfies Condition (*). Although generating DSR graphs for each 
X € X gives stronger results, application of the result is simplified if a single graph is 
generated. 

Example. As an illustration of Condition (*) consider the following pair of 
matrix-sets (the first consists of a single matrix) : 



A = 



-111 
1 -1 
1 -1 



B = 



-a b 
-c 
d 



(4.1) 



Assume that a,b,c,d > 0. The associated DSR graph Ga,b is shown in Figure | 
Although the es-cycles 5*1 — i?i — ^2 — i?2 and i?i — 5*1 — i?3 — 5*3 appear to have S-to-R 
intersection, the are not compatibly oriented, and in fact Condition (*) holds. 

4.4. Determinant expansions and structures in SR and DSR graphs 

Here we develop some methodology relating terms in determinant expansions and 
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Fig. 4.2. The DSR graph for the pair A,B in \4. j|) above. The graph satisfies Condition (*) 



objects in SR and DSR graphs. The most important notions are signed subterms 
in determinant expansions which are in one-to-one correspondence with signed term 
subgraphs in SR or DSR graphs. 

Below, A and B are sets of n x to matrices, Gj( and Gb are the associated SR 
graphs, and G = Gj^^b is the associated DSR graph. 7 = [71,72, 7fc] C {!,..., to} 
and 5 — [(5i, (^2, . . . , C {!,...,«} are fixed, ordered sets of equal size. 

Signed subentries in matrix-sets. Let A = {A{x)\x G X}, with Aij : X ^ M 
defined by Aij{x) — {A(x))ij, and Aij = {Aij(x) \ x G X}. Since Aij is a scalar func- 
tion on X, we can split it into Aij = A^j + A~j where Afj{x) = maxj^^j (a;), 0} and 

A~j{x) = mm{Xjix),0}. We cah Afj = {Atj{x) \x £ X} and ^r. = {A~j{x) \xeX} 
signed subentries in A with the usual rules for addition and multiplication. Corre- 
sponding to each of A'lj and A^j is an edge in G^. The same methodology can be 
applied to unsigned entries in B. 

Terms in determinant expansions: Given a permutation a of 7, define: 

^ l<5| _^ l<5| 

Ta^ P{a)Y[A5,c.,, T„ =P(a)]^65.a. . 

2=1 i=l 

For an interpretation of these expressions, the reader is referred back to Section [3?T1 
T a is a, term in the expansion of ^[(5|7], while is the corresponding term in the 
expansion of S[i5|7]. From now on, where an ordered pair of matrix-sets {A,B) is 
concerned, objects with left-pointing arrows above them are derived from the first 
matrix-set, and objects with right-pointing arrows above them are derived from the 
second. 

Signed subterms in determinant expansions. Given a permutation a of 7, 
consider the term (possibly zero) in the expansion of S[(5|7]. We can think of Tq, 
as a sum of signed subterms which can be enumerated in various ways. For example, 
for j = 1, . . . , 2l''l, define r^^^ to be the integer j — 1 written out as a binary string 
with \5\ digits. Then define the subentries — B^,^. if r^^^ has a zero in jth 

place, and T'f^^^ = B^^^^ if r(j) has a one in jth place. Many of these subentries may 
be zero if the corresponding entries are signed. Now Tq'' — P{a)Y[fli 'P^Pa- ^"^^ 
T a = X^jLi Ta\ The same methodology can be apphed to give signed subterms 

Tq ' in ^[(5|7]. From now on a "signed subterm" in the expansion of a determinant 
will mean a signed subterm not identically zero. 

If all entries Bs-a^ are nonzero and signed, then there is a unique signed subterm 
corresponding to a, i.e. T a — t\P for some j. Thus we can refer to t')P without 
further comment, to mean a signed subterm corresponding to permutation a. 
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Edges in the DSR graph are associated with signed subentries. We can 

associate with any nonzero signed subentry T^^^^, in .4(5|7) an S-to-R edge in G{5\^), 

which we term E^Pi^ . Similarly given some nonzero signed subentry T^'^^^ in B{5\^) 

we get an R-to-S edge in G{5\^), termed -E^'^^.- If for some a, j, T'f^^, and t'^^^, are 
both nonzero, then by construction they have the same sign, and hence correspond 
to the same (undirected) edge in G{5\^), which will be termed Eg'^,- 

Signed term subgraphs. We can associate with any signed subterm Ta in 
;6[5|7] an S-to-R edge-set in G{S\"f) which will be called Ea K Similarly for any signed 
subterm Tq^ in ^[5|7] there is an R-to-S edge-set in G{6\^), Ea\ Such edge-sets 
are a generalisation of term subgraphs introduced in [2] and will be called signed term 
subgraphs. Note that when we associate a signed term subgraph with a particular 
signed subterm we import this subgraph from the DSR graph including all directions 
and edge-labels. Further note that it only makes sense to talk about a signed term 
subgraph of an SR or DSR graph G, when G is square. Given a square DSR graph G, 
we will say that a signed term subgraph i? in G bisects a cycle G if it contains one half 
of the edges in G. This is only possible if it contains one member of a disconnecting 
partition of G. 

The above notions and notation carry over to SR graphs, except that in this case 
we omit the arrows indicating directionality. So given a matrix-set ^((5|7) and asso- 
ciated SR graph G_4(5|^), Ta is the subterm in v4[(5|7] corresponding to permutation a 
of 7, Ta^ is the jth signed subterm in Tq, and Ea^ is the corresponding signed term 
subgraph of G^(5|-y). 

Many of the results to follow rely on the fact that the union of two signed term 
subgraphs, one with S-to-R direction and one with R-to-S direction, results in a set 
of cycles in the DSR graph: 

Lemma 4.3. Consider two oppositely directed, signed term subgraphs Ea and E ^ 

— ^frl ^ — is) 

in some subgraph of a DSR graph. The union Ea U E , regarded as a subgraph, 
consists of a set of vertex-disjoint components, each of which is either an isolated edge 
or a genuine cycle. 

— ^ir) ^ — (s) 

Proof. Each vertex in Ea U E j^' either has incident on it two directed edges 

(one from Ea^ and one from E^^^), or one undirected edge. Firstly, if two vertices 

in Ea U E are connected by an undirected edge, then there are no other edges 
incident on either vertex. Secondly, no two formal cycles can intersect, as then there 
would be a vertex with three edges incident on it. Finally, all formal cycles are genuine: 
consider any formal cycle consisting of edges [ei, 62, . . . , e2r]. By the definition of a 
signed term subgraph, it is not possible for a and e(i mod 2r)+i to belong to the same 
signed term subgraph, and thus the edges must have alternating S-to-R and R-to-S 
direction. □ 

4.5. Hov^r the conditions in Lemmas 13.31 and 13.41 can fail 

Following the notation and terminology in Section [3. 11 consider two sets of square 
matrices of the same dimension A — {A{x) \ x G X}, B = {B{x) \ x e X}. We refer to 
{A, B) as a failed pair if there is some xq € X such that = -^(xq), and = B{xq) 
satisfy det(^o)det(i3o) < 0. {Aq,Bo) will be termed a failed instance of {A,B). 
With this terminology, some pair (5, V) fail the conditions in Lemmas 13.31 and 13.41 
iff there is some S C {1, . . . , n}, 7 C {1, . . . , m} such that (iS((5|7), {—V){j\S)) are a 
failed pair. 
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Consider a failed pair {A,B) with failed instance {Aq,Bq). There must be a 
nonempty set of signed subterms 7^ in the expansion of det(^) such that for any 
T G T^, Tdet(^o) C K>o, and similarly a nonempty set of signed subterms 7g in 
the expansion of det(S) such that for any T G Tg, Tdet(_Bo) C M>o. We will call 
the pair (7^,7b) failed subterms. 

5. Graph-theoretic results 

All results in this section build towards a proof of Theorem 14. II Two key prelim- 
inaries, Lemmas 15.11 and 15.21 below, are minor generalisations of results in [5]. The 
proofs are similar to those in [2 , and can be found in Appendix ID] 
Lemma 5.1. Consider two sets of k x k matrices A and B, and the associated DSR 
graph G = Ga.js- Consider signed subterms Ta'' in the expansion o/det(^) and T^^^ 



in the expansion ofdet{B). If all cycles in Ea' U E'^' are o-cycles then Ta'^T^^'^ C 

Lemma 15.11 tells us that a pair of oppositely directed, signed term subgraphs in 
a DSR graph derived from oppositely signed subterms must contain in their union 
some e-cycles. For example, define the following pair of matrix-sets: 



A 













" -bi 




* 




-a2 


















* 




as 






* 




b3 


* 


* 


a4 








* 






-bi 



(5.1) 



with ai,bi > 0, and * indicating entries of arbitrary magnitude and sign. Consider 
the corresponding DSR graph G_a,b- The entries in A define an R-to-S signed 
term subgraph in G^^g (Figure [5. lb ), and the entries bi in B define an S-to-R signed 
term subgraph in G^^g (Figure [5Tb ) . As the corresponding terms in det(^) and 
det(S) have opposite sign, the union of these term subgraphs contains an e-cycle 
fFigure 15.1b ). 



a) c 
R4 



i?2 



i?4 

1 

^4 



-R4 

i ■ i 

i?2 5^4 ^ R2 



Fig. 5.1. a) A signed term subgraph in the DSR graph G_i,B corresponding to a signed subterm 
in Aet{A) in h5. b) A signed term subgraph in G_^,g corresponding to a signed subterm in det(B) 
in h5. c) The union of these term subgraphs. As the terms are oppositely signed, the union of 
the corresponding term subgraphs contains an e-cycle. Edge-labels have been omitted. 



Remark. As the DSR graph does not contain information on whether A and 
B are independent. Lemma 15.11 is actually stronger than it appears. If all cycles 
in i^a'' U e'^p'^ are o-cycles, then in fact fixing any A £ A,B <E B and defining 

r„ = -P(a) nLi and T p = ^'(/5) IlLi -^ift , we have T o/T p > 0. The same 
holds true for each other DSR graph result: if a result about a pair of matrix-sets 
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{A,B) can be inferred from the DSR graph G.^ g, then it is true not only for each 
pair {A{x),B{x)) {x £ X), but also for each pair {A{xi), t3{x2)) {xi,X2 G X). 
Lemma 5.2. Consider any set of k x k matrices A with associated SR graph G. Let 
a and (3 be permutations 0/ {1, . . . , fc} such that Ta ^ and xjf^ are signed subterms in 
the determinant expansion of A. Assume that E^^ U E'^''' contains exactly one cycle 

C, and this cycle is an es-cycle. Then Ta'^ + T^^'^ = 0. 

Lemma 15.21 shows that having es-cycles in an SR-graph means that some terms 
in a determinant expansion sum to zero. The result is not affected by the possibility 
that some edge-labels in the graph may be 00, but obviously such labels cannot occur 
in es-cycles themselves. 

Lemma 5.3. Consider a set of square matrices A and the associated SR graph G^. 
Assume that Ga has an es-cycle C . Let £c be the set of all signed term subgraphs in 
Ga which bisect C, and Tc the corresponding signed subterms in det(y^). Then 

TdTc 



Proof. All edges in C correspond to signed subentries in A. If Ec is empty 
then we are done. Otherwise consider E £ £c with corresponding signed subterm 
T. Construct the new signed term subgraph E = {E\C) U {C\E) with corresponding 
signed subterm T. Clearly E £ £c- Now E U E contains a single e-cycle G which is 
an s-cycle, so by Lemma [5.21 T + T = 0. All signed subterms in Tc pair off in this 
way, so EtsTc T ^ 0. D 

Note that the above lemma does not imply that det(^) = 0, as not all signed 
term subgraphs in Ga necessarily bisect G. 

Lemma 5.4. Consider a set of square matrices A and the associated SR graph Ga- 
Let C be any set 0/ edge-disjoint es-cycles in G. Consider all signed subterms in the 
expansion of det(^) such that the corresponding signed term subgraphs bisect some 
es-cycle in C. These terms all sum to zero. 

Proof Let C — {G^^\ . . . ,0'^''^. Define to be the set of aU signed term 
subgraphs which bisect G*^*^, let 7^ be the set of corresponding signed subterms, and 
let T = U*L]^7^. By Lemma [531 all terms in % sum to zero. 

We will refer to a pair of signed term subgraphs E and E in £i such that [E U 
E)\{E n E) = G^*-' as an i-pair. Now consider some signed term subgraph E G £i CiSj 
for some i,j. Because the G^-'^ are all edge-disjoint, we have the following partition 
of 

E={En gW) u (£; n g(j')) u {e\{g^'^ u g^^'))) . 

Let E — (E\G'^^^) U (G^*''\-E) be the other member of the i-pair corresponding to E. 
We know that (En G^^)) C (E\G(*)), so 

E n g(^') = {E\c^% n g(^' = e n c^^^ 

which makes it clear that E £ £j. So if one member of an i-pair is in £j, then so is 
the other. Since Si CiSj consists of i-pairs, it follows that £i\£j consists of z-pairs. By 
induction, if « > 1, then £i\ U^^^ £j consists of «-pairs. So, by Lemma 15.31 

E ^ = 0. 
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We have the partition T = Ti U (7^\Ti) U (T3\(Ti U 7^)) U ■ • ■ , and hence 



E^-E E ^ = Eo 

TGT ^=l z=l 



□ 

Corollary 5.5. Consider a set of square matrices A and the associated SR graph 
Ga- Consider any set C of edge-disjoint es-cycles in Ga, md assume that each 
signed term subgraph in bisects some es-cycle from C. Then all matrices in A are 
singular. 

Proof. This follows immediately from the previous result. □ 

The matrix-set and corresponding SR graph in Figure [5T2l provide an illustration 
of the previous results. 



1 








2 





















1 








2 




Fig. 5.2. A matrix-set and the corresponding SR graph, 
nitude and sign (and may be unsigned), and thus wavy edges 
single edges or pairs of edges, and have unknown edge-labels, 
signed term subgraph in the DSR graph bisects the es-cycle Si 
are singular. 



The entries * are of unknown mag- 
in the SR graph may correspond to 
However, it is easy to see that every 
— Rl — S4 — R4 , and thus the matrices 



Lemma 5.6. Consider a pair (A, B) of sets of k x k matrices, and the associated DSR 
graph G = Ga.b- Assume that {A,B) are a failed pair with failed instance (^Qi-Bq) 
and failed subterms (7^,7g) corresponding to sets of signed term subgraphs {£j[,£js). 
Choose some T ^ E with corresponding .signed term subgraph £ £b- Let C be 
the set of all e- cycles in G each of which lies in E ^ U E for some E £ Then 

1. C contains an e-cycle which fails to be an s-cycle, or 

2. C contains two e-cycles which fail to be edge- disjoint. 

Proof. Tj( and Tg are nonempty (by the definition of a failed pair), and so, by 
Lemma [Ol C is not empty. Assume the result is false, i.e. C consists of edge-disjoint 
es-cycles G'''^\ . . . , G'^^\ All edges in all cycles in C correspond to edges with the 
same edge-labels in G^ (the SR graph corresponding to A): otherwise some edge 
must be imported only from Gg, and hence carry an edge-label of 00, causing some 
cycle in C to fail to be an s-cycle. Let Si be the set of all signed term subgraphs in 
Gjs, which bisect C^'\ with corresponding signed subterms %, and let T — UT^. By 
construction, T D X4- By Lemma |5.4[ the sum of terms in T is zero. Since ^0 is 
nonsingular, there must be some signed subterm T such that rdct(Ao) C M>o. 
By definition, T eTj^, contradicting the fact that T D 7^. □ 

The next result is about the geometry of subgraphs of DSR graphs constructed 
as the union of exactly three signed term subgraphs: 

Lemma 5.7. Consider a square DSR graph G containing .signed term subgraphs, Ei, 
E2 and E3. Assume that there is a cycle G in E i\J E 2 and another, distinct, cycle 
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D in E lU E3. C and D must either he edge and vertex-disjoint, or must have S-to-R 
intersection. _^ ^ 

Proof. Each edge in C lies in exactly one of Ei or E2, and similarly each edge in 
D lies in exactly one of Ei or E^. Consequently, any edge in CDD lies either in Ei 

> i 

or in n i?3, but not in both. Thus C and D each have a natural orientation and 
any edge in C n D has the same C-orientation and /^-orientation (i.e. C and D are 
compatibly oriented). Define G123 = EiUE2l-iE3, and ignore signs and labels on 
edges. No vertex from 6*123 can have more than three edges from G123 incident on it, 
so it is impossible for C and D to have a vertex but no edges in their intersection as 
then this vertex would have four edges incident on it. If C and D are vertex-disjoint 
then we are done. 

So assume that C and D share a vertex and choose any such vertex vq. Since 
C and D are distinct, following C backwards from vq we must come to a first vertex 
Vj (possibly vq) with three edges from G123 incident on it, two incoming edges (in 
C\D and D\C), and one outgoing edge (in C Ci D). Since E i, E2 and E^ are term 
subgraphs, the outgoing edge must lie in E'l, while the incoming edges must lie in 
E2 and i?3. Similarly following C forwards from vq we must come to a vertex Vk 
(possibly Vq) with three edges from G123 incident on it, one incoming edge (in CClD) 
and two outgoing edges (in C\D and D\C). This time the incoming edge must lie in 
El, while the outgoing edges must lie in E2 and E3. (see Figure [531) . Since it starts 
and ends with an edge from E 1 , the path in G n D from Vj to Vk has an odd number 
of edges. As vo was an arbitrary vertex in G n £*, G and D have S-to-R intersection. 




Fig. 5.3. A portion of the DSR graph corresponding to the situation in Theorem \ 5. 7| Because 
any vertex with three edges incident on it must have one from E\, one from E2 and one from E3, 
this forces any component of the intersection between C and D to be of odd length. 



As an immediate corollary of Lemma 15.71 
Corollary 5.8. Consider a DSR graph G with a signed, directed term subgraph Ei. 
Let E2, . . • , Ek (k > 2) be a set of oppositely directed term subgraphs in G. Define 
the set of e-cycles C in G as follows: a cycle C is in C iff C C {E i U E j) for some 
j £ {2, . . . , fc}. Then any two cycles in C are either edge and vertex-disjoint, or have 
S-to-R intersection. ^ 

Proof. Consider two cycles, G and D in C. Let Ei be some signed term subgraph 
such that EiUEi D G and Ej be some signed term subgraph such that EjUEi D D. 
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If i = J, then clearly C and D must be edge and vertex-disjoint since no vertex in 
Ei U El has more than two edges incident on it. Otherwise, apply Lemma 15.71 to 
El, Ei and E j to get that C and D must either be edge and vertex-disjoint or must 
have S-to-R intersection. □ 

As an illustration of Lemma [5771 and Corollary [5T8] consider the following matrix- 
sets (one of which happens to consist of a single matrix) : 





' 1 





1 " 




a 





" 







1 


1 


, B = 





b 







1 


1 


1 










c 



(5.2) 



Assume that a, 6, c > 0. A and B are a failed pair since det{A)det{B) = —abc. Define 
Ti = B11B22B33 = abc, T2 = -A11A23A32 ^-1, and_ T3 = -^13^22^31 = ^1, 
with corresponding signed term subgraphs E i, E2 and E 3. The subgraphs E lU E2, 
EiU E3 and E i (J E 2 U E 3 of the DSR graph are shown in Figure [3731 Each of 
El U E2, El U E3 contains an es-cycle. As these are not edge and vertex-disjoint, 
they have S-to-R intersection, consisting of the edge S3 — R3. 



Si 



EiU E2^ 



Ri 




Si <- 



EiU E2U E3 ^ 



Ri 



Si <- 



EiU E3 = 



■R3 



Ri 

1 



■i?2 



■R3 



■^S3 



R2 



Fig. -5.4. Subgraphs of the DSR graph for the pair A,B in i5.2{) above. By Lemma \ 5. 7] the 
distinct es-cycles are forced to have S-to-R intersection. 



We now come to the proof of Theorem 14.11 which states that if an interaction 
network defined by matrix-sets S and V fails the conditions in Lemma 13.31 then the 
associated DSR graph fails Condition (*). 

Proof of Theorem \4.1\ If S and V fail the conditions in Lemma l3.3l there are some 
nonempty S C {1, . . . , n} and 7 C {1, . . . , m} such that A = 5(^|7) and B = — V(7|^) 
are a failed pair, with nonempty sets of failed subterms (7^4, 7h). Recall that a term 
from and one from correspond to distinct signed term subgraphs in G((5|7). 
Choose some T ^ € Tjg with corresponding signed term subgraph i?* in G{S\j). Let 
£_4 be the signed term subgraphs corresponding to terms in T^. 

Let C be the set of e-cycles defined as follows: an e-cycle C is in C iff it lies in 
E^,\J E for some E G £j[. By Lemma [5.61 C must either contain an e-cycle which 
fails to be an s-cycle or two es-cycles which fail to be edge-disjoint. In the latter case, 
by Corollary [5711 these must have S-to-R intersection. In either case Condition (*) is 
failed. □ 
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6. Examples 



6.1. A system with three variables 

Consider the following dynamical system on M.^: 

±1 = /i(a;i,a;2) - gi(xi) 

i2 = f2ixi,X2,X3) - q2{x2) (6.1) 
^3 = f3{xi,X3) - qsixs) 

where > 0, and moreover |^ < 0, and 1^, 1^, |^ > 0. Without 

further information, there is no obvious decomposition of [/i, /2, fsf" except for trivial 
decompositions. Defining Vij = \dfi/dxj\ we have the Jacobian 





-I'll 


-V12 


" 




-Vll 


~Vi2 


" 




' 1 
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J = 


V21 


-V22 






V2I 


-V22 


V23 







1 







V31 





~V33 




V31 





~V33 










1 



Using this decomposition we can construct the DSR graph shown in Figure [Ol which 
contains no e-cycles. From this, we deduce that System \6A] is injective on M.^ (or any 
rectangular subset of K.^). This example shows that the trivial decomposition can be 
useful. The relationship between DSR graphs corresponding to trivial decompositions, 
and the usual interaction graphs [31 13] is discussed further in the concluding section 
and forthcoming work. 

xi Ri ^ X3 



i?2 ^2 < R3 

Fig. 6.1. The DSR graph corresponding to the trivial decomposition J = JI for System (6. 
All edge-labels are 00 and have been omitted. By inspection, the graph contains no e-cycles and 
hence satisfies Condition (*). As System 16. l\l has full outflows, it is thus injective. 



Remark. In this example, if ^j^, ||^, ||| < 0, then since the DSR graph contains 

no e-cycles at all, this means that even without outflows the system has p(~) Jacobian 
forbidding multiple equilibria. We do not pursue this here. 

6.2. A single very simple reaction 

Consider a single chemical reaction involving two substrates A ^ B. Denoting 
the concentration of A by a, that of B by 6, the reaction rate by v(a, 5), and ignoring 
any inflows and outflows, we get the dynamical system 



F{a,h) = 



-v{a, h) 
v{a, b) 



Defining S — [—1, 1]-^ to be the stoichiometric matrix of the system, gives the natural 
decomposition F{a, h) = Sv{a, h). However we could also choose to write F = F o id, 
where id(-) is the identity on R^. These two choices give us the following two product 
forms for the Jacobian DF{a, b): 



DF{a,b) 



[ 



Vb 



and DF{a, h) = 



-Va 


-Vb 




' 1 


" 


Va 


Vb 







1 
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where = ^ a-nd ^fc = ff- With the natural assumption that Wa > and vi, < 0, 
corresponding to the two products are the DSR graphs shown in Figure [ 



A 



Ri 



B 



A 



i?2 



Ri 



B 



Fig. 6.2. Left: The natural DSR graph for the reaction A B constructed from the decompo- 
sition F = Sv. The graph contains no cycles. Right: The DSR graph graph constructed from the 
decomposition F = F o id. The graph contains an e-cycle which fails to be an s-cycle. 



The natural DSR graph has no cycles and it is immediate that the system is 
injective provided both substrates are subject to degradation/outflow. On the other 
hand the DSR graph obtained from the decomposition F — F oid contains an e-cycle 
which is not an s-cycle, and cannot be used to draw such conclusions. This example 
illustrates the importance of the choice of decomposition, discussed further in the 
conclusions H 

6.3. Corollary [m 

Consider the result in Corollary 13 . 61 which gave submatrices and subgraph shown 
in Figure [6731 The subgraph has an e-cycle which fails to be an s-cycle. 

oo 

^ ^ Rk 

oo > 



" Sj, 






Sli Sik 




_ X Vki _ 



Fig. 6.3. Submatrices and subgraph of the DSR graph corresponding to the situation described 
in Corollarv \3.(A Wavy edges correspond to edges of unknown sign and label or possibly edge-pairs. 
The subgraph of the DSR graph has been drawn for Su = Vij = V^j = 0. Allowing these to take 
nonzero values simply makes some directed edges undirected. Clearly at least one of the cycles 
connecting Ri, Sj, Rk and Si must be an e-cycle which fails to be an s-cycle as it includes an edge 
with edge-label oo. 



The fact that the DSR graph contains an e-cycle that fails to be an s-cycle does 
not in itself prove that the system must fail the conditions in Lemmas 13.31 and 13.41 
- for this the matrix-theoretic formulation is necessary. It is an easy matter to find 
counterexamples [2] illustrating that Condition (*) is not necessary to ensure that a 
system is Pq This example does illustrate however that with certain assumptions, 
certain motifs in a DSR graph may be sufficient to ensure that the conditions are 
failed, a fact which will be explored in future work. 

■^It happens, in this case, that if for each fixed value of (a, b) we generate a DSR graph based 
on the decomposition F = F o id, then these all satisfy Condition (*) because the e-cycle is in each 
case an s-cycle. However it is often the case that for one decomposition the DSR graph satisfies 
Condition (*) globally, while for another it fails Condition {*), both pointwise and globally. 
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6.4. A famous interaction network: the "repressilator" 

A famous example of a synthetic oscillator is the repressilator described in |16j . 
The system is described by the differential equations: 



rhi 


= ao - 


^ f{p3) 


— mi 


7712 


= ao - 


^ fiPi) 


- m2 


TO3 


= ao - 


^ f{p2) 


- 7713 


Pi 






- Ppi 


P2 






- PP2 


P3 






- PP3 



where ao,/3 > and / is a decreasing function. The matrices S, — V"^ and the DSR 
graph G = Gg^^r for the system are shown in Figure [6^ With the given outflow 
conditions the system is injective on any rectangular domain as G consists of an o- 
cycle. Injectivity is also easily derived using the trivial decomposition, or using results 
from [31 S]. 
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Fig. 6.4. The matrices S and — V^, and the resulting DSR graph for the repressilator. The 
guantities va are all positive. The DSR graph consists of a single o-cycle. Edge-labels have been 
omitted. 



6.5. Reaction systems with unknown influences 

We present a reaction network involving two reactions and four substrates which 
is an example of a system which (with outflow) is injective for all kinetics despite the 
fact that we don't know whether an influence is activatory or inhibitory. Intuitively, 
this is because it has hierarchichal structure, decomposing into two subsystems, each 
well-behaved, and with only one-way influence between the two. The system, with 
corresponding matrices S and — V"^, and resulting DSR graph, is shown in Figure [^31 
As the graph contains no cycles, injectivity of the system with outflows is immediate. 

6.6. The importance of reversibility 

Consider the reaction systems and their respective DSR graphs shown in Fig- 
ure 16. 6[ drawn for any kinetics satisfying the NIC condition [2]. In each case A, a 
product of the flrst reaction, activates the second reaction. In Figure 16.6b . the flrst 
reaction is assumed to be reversible. The associated DSR graph has an e-cycle which 
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A + C ^B- 
C^D 



C 



-1 

1 

-1 

1 



-vii 

Vl2 

X -V23 

V24 




Fig. 6.5. A reaction system with matrices S and — V^, and the DSR graph. C modulates the 
first reaction with influence of unknown sign. The guantities Vij are all nonnegative (as expected for 
any reasonable kinetics), while x is unsigned. Edge-labels have been omitted. There are no cycles 
in the graph, and so the system is a Pg ' system. 



fails to be an s-cycle and so cannot be used to make injectivity claims. In Figure [6?6b 
the first reaction is assumed to be irreversible (in the strong sense that the products 
of the reaction cannot influence the reaction rate). In this case the DSR graph con- 
tains no cycles at all, and so the system with outflows is injective. Intuitively, this is 
because in the flrst case A and B may be able to activate their own production, while 
in the second they cannot. 

a) b) 

C^A+B C^A+B 

[ i 

B ^ D B ^ D 
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■Ri- 



A 



C 



A 



B 



■R2 



D 



B 



■i?2 



D 



Fig. 



Two reaction systems and their DSR graphs. All edges except for those labelled oo 

(-) 



have edge-label 1. // the first reaction is irreversible, then the system is a Pq system. 



6.7. Back to the TCA cycle 

We return to the DSR graphs of the TCA cycle model in [M] shown in Figure [TTTl 
A couple of points are noteworthy. First, several of the reactions appear as reversible 
in the DSR graph (in particular the production of fumarate and oxaloacetate) because 
although the functional forms given in [M] may only permit the reaction to proceed 
in one direction, the reaction rates are inhibited by the reaction products. Secondly, 
given a conserved pair such as NAD/NADH, one can write [NAD] = NADtot - [NADH] 
for some constant NADtot to eliminate NAD from the system and simplify the DSR 
graph. Subsequently all claims refer to the system with some fixed value of NADtot ) 
i.e. on some fixed stoichiometric class. This procedure works well with a conserved 
pair, but is not necessarily helpful with more general conserved quantities, possibly 
introducing a number of additional edges into the DSR graph. 

At the first stage of construction, shown in Figure II. lb . the DSR graph has a 
single es-cycle and fulfils Condition (*). Adding in the inhibition by oxaloacetate 
of succinate dehydrogenase (Figure II. lb ) introduces a single o-cycle with undefined 
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stoichiometry. It does not introduce an e-cycle with undefined stoichiometry, precisely 
because some of tlie reactions are irreversible (in the strong sense), indicating again 
the importance of considerations of reversibility which would be missed by analysis 
using the SR graph alone. Introducing NADH adds a single S-vertex and three extra 
positive edges (with edge-labels 1) to the DSR graph (Figure [Lib ). All these are 
undirected whether or not the reactions are reversible in the chemical sense, because 
an increase in the level of NADH is equivalent to a drop in the level of NAD, and so 
affects the reaction rate. This adds a great number of cycles to the graph, but they 
are all o-cycles. Finally, adding the AAT-catalysed reaction (Figure [TTTH ). leads to 
the creation of yet more cycles. Now the graph violates Condition (*), both because 
it contains e-cycles which fail to be s-cycles (e.g. OAA-i?6-FUM-i?7-MAL-i?8-NADH- 
i?3-aKG-i?9-OAA) and because it contains es-cycles with S-to-R intersection. 

Note that even at this stage some model quantities, in particular ADP, calcium, 
and membrane potential are omitted from the analysis. This example illustrates that 
even in models where we cannot positively rule out multiple equilibria, by identifying 
cycles which cause Condition (*) to be violated, we can speculate on the mechanisms 
by which multistationarity may arise. 

7. Discussion and Conclusions 

We have shown that the DSR graph, an object closely related to the pathway 
diagrams and interaction diagrams drawn by applied scientists in various fields, can 
be used to rule out the possibility of multiple equilibria or multiple nondegenerate 
equilibria. Given a systems written in one of the forms 

X ^ f{vix)) ~ Q{x) or X ^ f{v{x)) - Qeix) , 

we construct the DSR graph for the associated systems x — f{v{x)), and check 
whether Condition (*) holds for this graph. If it does then we know that all of 
f{v{x)) — Q{x) are injective and that all of i = f{v{x)) — Qg{x) forbid multiple 
nondegenerate equilibria on the relative interior of any invariant affine set. 

Algorithmic development is clearly important at this stage. Simple, widely avail- 
able, computational tools to test whether Condition (*) holds for a DSR graph will 
open up the possibility of routine preliminary model analysis: modellers working with 
qualitative models will be able to check whether their models, as a result of struc- 
ture alone, forbid multiple equilibria, before parametrisation and simulation. We are 
currently working on creating user-friendly software that performs this test. 

It is worth remembering that matrix-theoretic approaches give sharper results 
than the graph-theoretic ones. As mentioned in the text. Lemmas 13.31 and 13 . 41 simplifv 
in a number of special cases where we have additional knowledge of a system (e.g. 
reversibility of all processes). This specialisation is a task for future work. 

A number of special cases can be treated via minor variants on this analysis. 
Some systems on the nonnegative orthant, for example population models, generally 
have boundary equilibria, and one question is whether they admit multiple positive 
equilibria. The general theory developed here can be applied in this context, for 
example to a Lotka-Volterra type model for the evolution of n species of the form 
Xi = XiGiix) [22[ I23j where x — [xi,. .. ,Xn]'^ is the nonnegative vector of species 
concentrations, and the functions Gi{x) encode information on the interactions be- 
tween species. Positive equilibria must satisfy Gi{x) — 0, and hence we are interested 
in injectivity of G{x) — [Gi{x)^ . . . ,G'„(x)]^ in the interior of the positive orthant, 
which can be treated using the techniques developed here. 
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The example in Section 16.11 illustrated that even trivial decompositions can be 
useful. In fact, in forthcoming work, we will show how given an arbitrary dynamical 
system x = F{x) on a rectangular domain, and using the trivial decomposition F = 
F o id, it is possible to obtain stronger results than the well-known interaction graph 
results on injectivity (results in |4j for example). However, the example in Section [6. 21 
illustrated the important and subtle point that in some cases the strongest results are 
obtained by choosing the "correct" decomposition of the functions we hope to show are 
injective. This is true for both the matrix-theoretic and the graph-theoretic results. 
The following difficult question remains open: given an arbitrary set of functions, how 
do we choose the decomposition which will give us the sharpest results on injectivity 
(assuming that such a choice exists)? 

Finally, the importance of cycle structure in SR and DSR graphs goes beyond 
questions of injectivity. Cycle structure in SR graphs is, for example, linked closely 
to monotonicity [9] (several results in also have immediate graph-theoretic in- 
terpretations). The close relationship between cycle structure and the possibility of 
complex behaviour has been understood for some time for models whose structure 
can be represented via an interaction graph [5] , but the corresponding results for gen- 
eral interaction networks represented via SR/DSR graphs (or possibly other variants 
of these ideas) are less complete. Such study has important implications for model 
caricature and model simplification. 
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Appendix A. Definitions and notation. 

In each of the the definitions below, M is an n x m matrix, S C {1, . . . , n}, and 
7 C {1, . . . , m} with \S\ — \ 

M{d\^) is the submatrix of M with rows indexed by S and columns indexed by 
7. A principal submatrix of M is a submatrix of the form M{S\5). Determinants 
of submatrices of M are termed minors of M. Principal minors are determinants 
of principal submatrices. 

P matrices are square matrices all of whose principal minors are positive. They 
are by definition nonsingular. If —M is a P matrix, then M is a p(^) matrix. Each 
k X k principal minor of a p(~) matrix has sign (—1)*^. 

Pq matrices (so termed in [24]) are matrices in the closure of the set of P 
matrices. These are matrices all of whose principal minors are nonnegative. We will 
term M a Pq ^ matrix if —M is a Pq matrix. A k x k minor of a Pq matrix is 
either zero or has sign (—1)*^. The zero matrix is a Pq and Pq matrix. 

Qualitative classes. A matrix M determines the qualitative class Q{M) [2S] 
consisting of all matrices with the same sign pattern as M. Explicitly, Q{M) consists 
of all matrices X with the same dimensions as M, and satisfying Mij > ^ Xij > 0, 
Mij <0 ^ Xij <0 and Mij = ^ Xij = 0. 

Qo{M) is the closure of Q{M). Explicitly, Qo{M) consists of all matrices X with 
the same dimensions as M such that Mij > ^ Xij > 0, Mij < => Xij < and 
Mij = ^ Xij = 0. 

A square matrix M is sign nonsingular if all matrices in Q{M) are non- 
singular. A matrix-set A4 has signed determinant if either det(A^) = 0, or 
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det{M) C (-00,0) or det{M) C (0,oo). 



Appendix B. Multiple nondegenerate equilibria in systems with con- 
served quantities. 

The discussion in this appendix is closely related both to Theorem 2 in pS' and 
to Proposition 1 in 4 . 

Consider a dynamical system x = fix) defined on X C M". An equilibrium p of 
the system is nondegenerate if Df{p) (the Jacobian at p) has no zero eigenvalues. 
When X is foliated by invariant sets we are generally interested in behaviour on one of 
these invariant sets. Assume that i — f{x) preserves some function E : R" — *■ 
(1 < Z < n), so that E{x) = along trajectories. Since any level set of the form 
Ec = {x\E(x) = C} in invariant, the pertinent question is whether the system 
restricted to EcC\X admits multiple nondegenerate equilibria, i.e. whether the system 
X = f{x) can have two equilibria p,q ^ [Ec H X) such that Df{p) and Df{q) have 
no eigenvectors corresponding to zero eigenvalues, tangent to Ec- 

Assume that A is a rectangular subset of R" . We already know that if (|AoD is a 



Pq •* system, then ( A+ ) is injective on X. We now show that if ( A+ I is injective on 



X , then dAel ) is incapable of multiple nondegenerate equilibria in the following sense: 
let Ec be any invariant afRne subspace of K". Then the system restricted to Ec n X 
can have no multiple nondegenerate equilibria in ri^Ec H X), the relative interior of 

Ecnx. 

In order to show this, we need the following basic persistence property of nonde- 
generate equilibria: if a vector field / on some smooth manifold S (possibly with 
boundary) has a nondegenerate equilibrium p 6 int(5) then any vector field g on S, 
close to / in the topology, has an equilibrium p , close to p, in int(iS). This follows, 
for example, by isolating p in a sufficiently small closed neighbourhood Up C int(tS) 
and using the invariance of the Brouwer degree in int([/p) under small perturbations 
of/. 

Lemma B.l. Consider a dynamical system x = f{x) on X C R", where X is a 
convex, forward invariant set. Assume that 

1. Given any A > 0, all functions of the form f{x) — \x are injective on X . 

2. There is an affine subspace Ec C M" such that Ec H X is invariant for 
X = fix). 

Then the system x — fix) restricted to Ec H X cannot have more than one nonde- 
generate equilibrium in liiEc HA). 

Proof. The result is trivial if EcClX is empty or if Ec is a point. Otherwise, either 
Ec = R", in which case Ec^X = A, or Ec has dimension n~l for some 1 < / < n — 1, 
and there is a surjective linear function E : M" — > R', and a vector C G R' such that 
Ec = {x I Ex = C}, and for x £ Ec, Efix) = 0. In cither case, choose any vector 
k e Ec r\ X. Define the systems x = fix) + A(fc - x) = Fx{x) with A > 0. By 
convexity of A and the fact that A is forward invariant for x = fix), we have that 
A is forward invariant for x — Fxix). Further, x — F\ix) leaves Ec H A invariant: 
this is obvious when Ec = R"; when Ec + R", we have Efix) + \iEk ~ Ex) = Q 
when X G Ec- By assumption fix) — Xx, and hence Fxix), are injective on A. Now 
suppose that x = fix) contains two nondegenerate equilibria in viiEc HA). Since, 
for small A, F\ix) is close to fix), by the above persistence arguments, the system 
X = F\ix) must contain two (nondegenerate) equilibria on liiEc H A) contradicting 
the fact that Fx are injective. □ 

Remark. In Lemma lB.ll rather than choosing equilibria, i.e. values of x satisfy- 
ing fix) — 0, we can consider values of x satisfying fix) = c for any c. Theorem 4w 
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in [50] implies that if X is an open rectangular subset of K", and, for some 0, \Ng\ is 
a Pq ' system with nonsingular Jacobian everywhere in X, then it is injective on the 



interior of X. This result can be seen as a special case of Lemma [B. II 

We remark that extensions of Lemma IB.ll to allow nonlinear integrals E are also 
possible - the key question is whether given some fixed level set Ec we can construct 
a family of systems x — F\{x) which preserve Ec, are injective for A > 0, and 
such that Fi^{x) — f{x). We do not pursue this here. 

Appendix C. Further discussion of DSR graphs. 

DSR graphs can be understood via possible motifs. Consider some species with 
concentration 5* which participates in some interaction with rate R. In Figure IC.ll 
all six possible single-edge connections between an S-vertex and an R-vertex are illus- 
trated. There are nine possible double edge connections between two vertices: These 
are created by taking any of the motifs in Figure IC.lb . b, or c, and combining with 
any motif from Figure [C. Id . e, or f. Two common ones are shown in Figure [C. 21 and 
their meanings are described in the caption. 

a) h) c) 

S > R S < R S R 

d) e) f) 

S > R S « R S R 

Fig. C.l. All single edge connections between an S-vertex and an R-vertex in a DSR graph. 

a) Increases in S decrease R (S inhibits the interaction), but S is unaffected by the interaction. 

b) The interaction increases S, but its rate is not affected by S. c) Increases in S decrease R and 
the interaction increases S. d) Increases in S increase R (S activates the interaction), but S is 
unaffected by the interaction, e) The interaction decreases S, but its rate is not affected by S. f) 
Increases in S increase R and the interaction decreases S. 



a) ^ ^ b) ^ ^ 

^ R ^ - ^ , R 

Fig. C.2. Two common double edge connections between vertices in a DSR graph, a) Increases 
in S may cause either an increase or a decrease in R, but S itself is unaffected by the interaction, b) 
Increases in S may cause either an increase or a decrease in R, and S is decreased by the interaction. 



The DSR graph is a natural amalgamation of directed versions of two SR graphs 
Gs and G_yT, associated with two sets of matrices, S and V. All the matrix- 
theoretic results treat S and V symmetrically, and ideally this would be reflected in 
the construction of the DSR graph. However, this would involve introducing two sets 
of edge-labels onto the graph. In our construction, for simplicity, only one set of 
edge-labels is imported from Gs, based on the practical fact that we often expect S, 
rather than V, to have constant entries and hence nontrivial edge-labels. However, it 
is straightforward to extend the treatment, allowing two sets of edge-labels, if required 
in applications. 

For a CRN, the DSR graph encodes information about irreversibility of reactions, 
and modulation by quantities which do not formally participate in the reaction. For- 
mally, it is the usual SR graph with the following modifications: some edges have 
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become directed, some short cycles have been replaced with single directed edges 
of defined sign, but edge-label oo, and some edge-labels on short cycles have been 
changed. If an o-cycle (resp. e-cycle, resp. s-cycle) survives, then it remains an o- 
cycle (resp. e-cycle, resp. s-cycle). It is clear that results for CRNs obtained using 
DSR graphs are sharper than those using SR graphs. 

Appendix D. Generalisations of results from [2]. 

The results here are self-contained, but [2j contains a more detailed discussion of 
the relationship between permutations of ordered sets and cycles in SR graphs. 
Lemma D.l. Consider a permutation a of some ordered set. Let a be written as the 
product of disjoint cycles from some set C. Let 6 — Ucgcc- Then P{a), the parity of 
a, is given by 

P(a) = (-l)|e|-|ci. 

Proof. This follows by writing any permutation as the product of disjoint cycles 
and noting that a fc-cycle is an even permutation if k is odd and vice versa. □ 

The following two lemmas are key results on the relationship between the signs of 
two signed subterms and the cycle structure of the union of the corresponding signed 
term subgraphs, first for SR graphs, and then for DSR graphs. The proof of the first 
result is presented in full for completeness. 

Lemma D.2. Consider any set of k x k matrices A with corresponding SR graph Gj^. 
Consider any two signed subterms T^^ and T^*^ in the determinant expansion of A, 
corresponding to permutations a and (3 of {1, . . . , /c} and signed term subgraphs Ea^ 
and e'p^ in Ga- Then 

Ti'-^rf;^ C (-1)"^='M>o (D.l) 

(r) (s) 

where \Ce\ is the number of e-cycles in Ea U Ep . 
Proof. By definition 

rMT(^) = p(a)p(/3)nr(:)r/;). 

i=l 

Define 

k 

Z ^ P(a)P(/3) nsign(i?l;!)sign(£;g) . 

i=l 

Proving the lemma is equivalent to proving that Z = (—1)''''='. 

Let 9 be the set of indices for which ElJ, and Eli^\ are distinct edges in G^. When 
z G {1, . . . , k}\e, sign(i?(^|)sign(i?,g) = L So 

Z = P{a)P{p) nsign(£;l:!)sign(ii;g). 

Now the edge set 
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(r) ( s) 

consists precisely of the set of cycles in Ea U E^^ ' . These cycles are vertex-disjoint 
since no more than two edges from this set can be incident on any vertex. Let the set 
of o-cycles in this set be Co and the set of e-cycles be Ce, with C ^ Co^Ce- Associate 

(r) (s\ 

with each cycle c a Co^ Ce the corresponding index set c, i.e., i ^ c E^J, , E^^' G c. 
Thus corresponding to the sets Co and Ce are the sets of index sets Co and Cg. Since 
any two cycles are edge-disjoint, Co U Cg is a partition of 9, and we can define 

0o = U c, e,= U s with \9o\ = J2 1^1' i^-i = E 1^1 ■ 

Clearly 6 — 6oLI 6^- We can write 




So 



VceCo / VceCe / 



= P(a)P(/3)(-l)l«°l+l''=l-l'^°l 
= P(a)P(/3)(-l)l'^l-l'^°l. 

Applying Lemma ID. II to l3 o gives us that 

P(a)P(/3)=P(/3oa-i)^(-lf-^ 

where g is the number of cycles in (3oa~^ and p is half the number of elements in these 
cycles. Now there is a one-to-one correspondence between nontrivial cycles in /5 o 

(r) (s) _1 

and long cycles in Ea U E^ , with a cycle of length I in f3 o a corresponding to a 

(r) (s) ' fr) is) 

cycle of length 21 in i?Q U E'^ . Let q be the number of short cycles in Ea U Ejj , 
so that q + q — \C\. Since there are precisely two edges in a short cycle, q is also half 
the number of edges in short cycles in Ea U Ep , so that p + q = \0\. This gives us 
that p — q=\6\ — |C|, so that: 

P(a)P(/3) = (-l)l^l-ici . 

Completing the argument: 

Z = (-l)l«l-|ci(_i)|e|-|Co| ^ (_j^)2|e|-|c|-|c„| ^ („;l)I'^I+I'^°I = . 

This proves the result. □ 

In the above result, it is perfectly possible to have a — (3 and/or r — s. Further 
the result is independent of edge-labels, so some of these may be oo. 
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For DSR graphs this resuh becomes: 
Lemma D.3. Consider two sets of k x k matrices A and B, and the associated DSR 
graph G. Consider any two signed subterms Ta^ in det(^) and T^^^^ in det{B), 
corresponding to permutations a and (3 of {1, . . . , k} and corresponding to oppositely 
directed, signed term subgraphs Ea^ and E^^^ in G. Then 

T^T^^^ C (-1)I'^^Im>o (D.2) 

where \Ce\ is the number of e-cycles in Ea U E ^ . 

Proof. Note first of all that by Lemma [4.31 all formal cycles in Ea U Ej^' are 
genuine cycles and are vertex-disjoint. From here the proof proceeds identically to 
the result for SR graphs. □ 

We now prove that if two signed term subgraphs have only o-cycles in their 
intersection, then the corresponding signed subterms have the same sign. First for 
SR graphs: 

Lemma D.4. Consider any set of n x n matrices A. Let T^^^ and Tg'^ be signed 
subterms in the expansion o/det(^). If all cycles in Ea^ U Ep''^ are o-cycles then 

T^C^)t^(s) r- TIP 
la -L p C_ K>o. 

Proof. If all cycles in Ea ^LlE'^'-' are o-cycles, then \Ce\ — 0, and applying Eq. (jP.ip 
immediately gives T^'^p'f^ C K>o. □ 



The result for DSR graphs was stated as Lemma [Ol 



Proof of Lemma \5.1[ If all cycles in Ea'^ U E^g^ are o-cycles, then \Ce\ — 0, and 



applying Eq. ^^^2^ immediately gives Ti'^^T^'^ C M>o. □ 

Finally, we prove Lemma 15.21 which stated that when two signed term subgraphs 
of an SR graph contain in their union a single cycle, and this is an es-cycle, then the 
corresponding terms sum to zero. 

Proof of Lemma \5.'A Note that the fact that C is an es-cycle implies that none 
of the edges in C have edge- label oo, even if some of the other edges in G may have 
edge-label oo. Thus C is a long cycle. By definition 



As usual, let 6 be the set of indices for which E^^, and are distinct edges in G, so 
that {£'j-^^.}ige and are precisely the edges in the unique es-cycle C. Defining 



Ci = {El^^,}i(ze and C2 — {sf^'lige gives us a disconnecting partition of C. Since C 

an s-cycle, val(Ci) and val(C2) are defined and equal. Define Z = Hieli k}\e -^i 
We can write 



is 

(r) 



\ iee iee / 

i\[Tl:]+P{Poa-')\{T^A 

\iee lee ) 

P{a)Z (sign(Ci)val(Ci) + P{(3 o a-i)sign(C2)val(C2)) 
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l3 o a ^ can be written as a single cycle of length \0\, and so from Lemma ID.li 
P(/3oa-i) = I.e., 

Ti'') + T^^) = P(a)Z (sign(Ci)val(Ci) + (-l)l«l-isign(C2)val(C2)) . 
Since C is an e-cycle: 

sign(C2)/sign(Ci) = sign(Ci)sign(C2) = sign(C) = (-1)1^1. 

(r) ( s) 

Substituting into the expression for Ta + gives: 

T^O + T^^) = P(a) Zsign(Ci) (val(Ci) - val(C2)) . 
Since C is an s-cycle, val(Ci) - val(C2) = 0, giving T^*^^ + ri"^ = 0. □ 
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